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We extend the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model by introducing an effective four- 
quark vertex depending on Polyakov loop. The effective vertex generates entanglement interactions between 
Polyakov loop and chiral condensate. The new model is consistent with lattice QCD data at imaginary quark- 
number chemical potential and real and imaginary isospin chemical potentials, particularly on strong correlation 
between the chiral and deconfinement transitions and also on the quark-mass dependence of the order of the 
Roberge-Weiss endpoint. We investigate the influence of the entanglement interactions on the location of the 
tricritical point at real isospin chemical potential and on the location of the critical endpoint at real quark-number 
chemical potential. 
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I. INTRODUCTION 

An important query on the thermodynamics of quan- 
tum chromodynamics (QCD) is whether the chiral-symmetry 
restoration and the confinement-to-deconfinement transition 
take place simultaneously or not. If the two transitions do 
not coincide, phases such as the constituent quark phase HHH] 
or the quarkyonic phase [3,4] may appear. 

If the chiral and deconfinement transitions are of first or- 
der, discontinuities appear simultaneously in their (approxi- 
mate) order parameters, that is, the chiral condensate a and 
the Polyakov loop <P |5l |g]. Furthermore, if a nontrivial 
critical endpoint (CEP) exists at finite temperature (T) and 
quark-number chemical potential /i q [7], susceptibilities of a, 
<P and other quantities diverge simultaneously |8|]. This in- 
dicates a coincidence of second-order phase transitions. At 
zero fi q , the chiral and deconfinement transitions are found 
to be crossover Hence, there is no a priori reason why 
the two transitions coincide exactly. Actually, in lattice QCD 
(LQCD) simulations at zero chemical potential B9l 4l IQ , there 
is a debate as to whether the transitions really coincide or not; 
see Ref. 11211 and references therein. LQCD simulations are 
far from perfection at real /i q because of the well-known sign 
problem 1131] ■ F ortunately, LQCD data are available at imagi- 
nary /i q [ 14-22] and real and imaginary isospin chemical po- 
tential fj,\ 123142511 . since LQCD has no sign problem there. 
The data show that chiral and deconfinement crossover tran- 
sitions coincide within the numerical accuracy. Since there 
is no general reason for exact coincidence between crossover 
transitions, it is natural to think that the chiral and deconfine- 
ment crossover transitions nearly coincide as a result of strong 
correlation (entanglement) between a and <P. We investigate 
this possibility in the present paper. 

As an approach complementary to first-principle LQCD, 
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we can consider effective models such as the Nambu-Jona- 
Lasinio (NJL) model 0, H and the Polyakov-loop ex- 
tended Nambu-Jona-Lasinio (PNJL) model SEH-IS. The 
NJL model can describe chiral symmetry breaking, but not the 
confinement mechanism. The PNJL model is designed lO 
to make it possible to treat both the mechanisms. The PNJL 
model can reproduce results of LQCD at zero and imaginary 
£t q i47l - l50ll where LQCD has no sign problem. 
At imaginary /.i q = i9 q T, the grand canonical partition func- 
tion ZQc(O q ) of QCD is related to the thermodynamic poten- 
tial i?QCD as J?QCD(0 q ) = -Tln(ZGc(0q)), where q is a 
real parameter. Roberge and Weiss (RW) found l55ll that QCD 
has a periodicity J?QCD(#q) — ^QCD($q + 2nk/3), showing 
that i7qcD(^q + 27rfc/3) is transformed into fiQCD^q) by the 
Z3 transformation with integer k. This means that QCD is in- 
variant under a combination of the Z3 transformation and a 
parameter transformation 9 q — » 6 q + 2/c7r/3 ll47il48ll . 

q -> Uq, A, -> UA V U- X - i/g{d u U)U~\ 
q -> 6» q + 27rfc/3, (1) 

where U(x,t) are elements of SU(3) with U(x,/3 = 1 /T) = 
exp(— 2iirk/3)U(x, 0), q is the quark field and A v is the 
gauge field. We call this combination extended Z3 transfor- 
mation. Thus, ^QCD^q) has the extended Z3 symmetry, and 
hence quantities invariant under the extended Z 3 transforma- 
tion have RW periodicity 11471 l48ll . At the present stage, the 
PNJL model is only a realistic effective model that possesses 
both extended Z 3 symmetry and chiral symmetry lt47l l48ll . 
Furthermore, the PNJL model can reproduce the first-order 
RW transition JH] that occurs at q = (2k + l)r/3 when T is 
larger than some critical temperature Tg. This property makes 
it possible to compare PNJL results with LQCD data quantita- 
tively at imaginary /i q . A current topic at imaginary /i q is what 
the order of the RW transition is at the endpoint T = Tg. The 
recent LQCD simulations show that it is first-order for small 
and larger quark masses, but the order is weakened and could 
be second order at intermediate masses lfl9ll20ll . 

In the PNJL model, the correlation between a and ^ is 
weak, so that the chiral and deconfinement crossover tran- 
sitions do not coincide without any fine-tuning of parame- 
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ters 14911 . For zero chemical potential, the scalar- type eight- 
quark interaction is necessary to obtain a coincidence be- 
tween the two transitions, and for imaginary /i q the vector- 
type four-quark interaction is needed 114911 . This fact indi- 
cates that a true correlation between a and is stronger 
than that in the PNJL model appearing through the covariant 
derivative between quark and gauge fields. Actually, recent 
analyses 11561 15711 ba sed on the exact renormalization-group 
(ERG) equation l58ll indicate that entanglement interactions 
between a and appear in addition to the original entangle- 
ment through the covariant derivative. 

In this paper, we extend the PNJL model by introducing an 
effective four-quark vertex depending phenomenologically on 
0. The effective vertex generates entanglement interactions 
between a and <P. The functional form of the entanglement 
vertex is determined by respecting the extended Z3 symmetry 
and the chiral symmetry. The strength of the vertex is deter- 
mined from LQCD data at imaginary /Lt q , and the validity of 
the model setting is confirmed for real and imaginary values 
of hi by comparing the model results with LQCD data. The 
new model is consistent with all LQCD data at imaginary ^ q 
and real and imaginary fij. Particularly, the new model can 
reproduce two phenomena simultaneously; one is the strong 
correlation between the deconfinement and chiral transitions 
and the other is the quark-mass dependence of the RW end- 
point predicted by LQCD very recently fl9l I2O1I . We also 
analyze the influence of the entanglement interactions on the 
location of the tricritical point (TCP) in fii-T plane and the 
location of the critical endpoint (CEP) in fi q -T plane. The 
present phenomenological approach is complementary to the 
ERG approach mentioned above. 

In Sec. II, we explain the PNJL model briefly and intro- 
duce an effective four-quark vertex depending on 0. In Sec. 
Ill, the new model with the effective vertex is applied to the 
imaginary /i q region and the real and imaginary /ii regions and 
compared with LQCD there. Sec. IV is devoted to a summary. 



II. PNJL MODEL 

We start with the standard two-flavor PNJL Lagrangian IL33L 

a 

£ =q{il v D u - m )q 

+ G s [(qq) 2 + [qi^Tqf] -U($[A],$[A}* ,T), (2) 



with 



r r p 

L(x)=Vexp i / drAi(x, r) 
L Jo 



(4) 



where V is the path ordering and A4 = iAq. In the chi- 
ral limit (mo = 0), the Lagrangian density has the exact 
SU(N f ) L x SU(N f ) R x J7(l) v x SU{3) C symmetry. The 
temporal component of the gauge field is diagonal in flavor 
space, because color and flavor spaces are completely sepa- 
rated in the present case. In the Polyakov gauge, L can be 
written in a diagonal form in color space 03311 : 

L = e # = diag(e i/3 ^,e Wb ,e w '=), (5) 

where <f> a = + <j) S /V3, (j) b = -fa + <t>s/V^ and <\> c = 
— {4> a + 0h) = —^4>s/V^>- The Polyakov loop is an exact 
order parameter of spontaneous Z3 symmetry breaking in pure 
gauge theory. Although Z 3 symmetry is not an exact one in 
the system with dynamical quarks, it still seems to be a good 
indicator of the deconfinement phase transition. Therefore, 
we use •P to define the deconfinement phase transition. 

Making the mean field approximation and performing the 
path integral over the quark field, one can obtain the thermo- 
dynamic potential J? (per volume), 



Q 



2N f 



d 3 P r 



(2tt) 3 



3E(p) 



+ i In [1 + 3(0 + ^-^(p)^-^^ + e' 3 ^'^} 
+ i In [1 + 3(<T + <£ e ^ + (p)) e -/^ + (p) + e -3^ + (p)] 
+ U M +U, (6) 

where a — (qq), S s = —2G s a, M = m + S s , Um = 
G s a\ E(p) = VP 2 + M 2 and E ± (p) = E(p) ± fi q = 
E(p) ± i8 q //3. On the right-hand side of ©, only the first 
term diverges. It is then regularized by the three-dimensional 
momentum cutoff A ll33l - l35ll . We use U of Ref. [38], which 
is fitted to a LQCD simulation in pure gauge theory at finite 
T 



U = T 4 



a(T) 



b{T) ln(l - 600* + 4(0 3 + 0* 3 ) - 3(00* 



(7) 



where q denotes the two-flavor quark field, mo denotes the 
current quark mass, and D v = d v + iA v — ifiqSg. Field A v 
is defined as A v = S%gA° a ^-, with gauge fields A v a , the Gell- 
Mann matrix A a , and the gauge coupling g. In the NJL sec- 
tor, G s denotes the coupling constant of the scalar-type four- 
quark interaction. The Polyakov potential U, defined in (0, 
is a function of Polyakov loop and its Hermitian conjugate 

0*, 



1 

= — TrL, 



0* = -TVZt, 



(3) 



with 



a(T) = ao + ax(^Y) 



(8) 



where the parameters are summarized in Table U In pure 
gauge theory, the Polyakov potential yields a first-order de- 
confinement phase transition at T — Tq. It is determined 
from pure gauge LQCD that T = 270 MeV. In QCD with 
two-flavor dynamical quarks at n q — 0, the PNJL model 
with the original value of Tq shows that the pseudocritical 
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temperatures of chiral and deconfinement crossover transi- 
tions are T a w 230 MeV and T$ w 215 MeV, respec- 
tively S, while full LQCD simulations |1 ED, SI show 
that T a ~ Ttp ~ 173 i 8 MeV. It follows from these re- 
sults that the relative difference A = \T a — T<p\/T a is about 
6% for the PNJL model and at most 10% for LQCD. Thus, 
for A the PNJL result is consistent with the LQCD data, but 
for the absolute values of T a and T$ the PNJL result is larger 
than the LQCD data. Therefore, we rescale T to 212 MeV 
in the PNJL model to obtain 7$ = 173 MeV. However, the 
PNJL calculation with the four-quark interaction only shows 
that T a = 216 MeV and A « 20% JH. The results on T a 
and A are not consistent with the LQCD results. This indi- 
cates that the entanglement between the chiral and deconfine- 
ment transitions is weak in the PNJL model. This problem 
will be discussed later and solved in this paper. 



do 


ffli 


a 2 


63 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov sector used 
in Ref. 1 3 8] . All parameters are dimensionless. 



The variables X = <P, <P* and a satisfy the stationary con- 
ditions 



dn/dX = 0. 



(9) 



The solutions of the stationary conditions do not necessarily 
yield a global minimum Q. There is a possibility that they 
yield a local minimum or even a maximum. We then checked 
that the solutions yield a global minimum when the solutions 
X(9 q ) are inserted back into ©. 

The thermodynamic potential fl of (O is not invariant under 
the Z3 transformation, 



<P(0 q ) -> <P(9 q )e 



;27rfe/3 



(10) 



although W of is invariant. Instead of Z3 symmetry, 
however, f? is invariant under the extended Z3 transforma- 
tion m, 



(11) 



This is easily understood as follows. It is convenient to in- 
troduce the modified Polyakov loop = e and = 
e~ i<P* invariant under the transformation ( TTTb . The ex- 
tended Z 3 transformation is then rewritten as 

<F(0 q )* -><F(0 q )% (12) 



and J7 is rewritten as 

n = 2Nf I Wy i 3E{p) + \ ln [1 + Me ~ mp) 



\ ln [1 + 3<P* e -^ p > + 3^e- 2 ^( p 'e- 



/3 



e -3PB(p) e - / 9 ( u B l 



(13) 



where /3/xb = 3/3/i q 



3i6» Q 



Obviously, j? is invariant un- 
der the extended Z3 transformation (fT2l . since it is a func- 
tion of only extended Z 3 invariant quantities, e 3iSq and X(= 
if - , a). The explicit 9 q dependence appears only through a 
factor e 3 ' q in ( TT3T >. Hence, the stationary conditions (O show 
that X = X(e 3te<1 ). Inserting the solutions back to ( fT3l . one 
can see that fl = 0(e 3ie «). Thus, X and Q have the RW 
periodicity, 

x(d q + ^) = x(e q ), and r?(^ q + ^) = r?(e q ), 

(14) 

while the Polyakov loop ^ and its Hermitian conjugate 
have the properties 



27Tfc N 

*{e q + — 

w¥r 



a i2nk/3 



(15) 



The RW periodicity is a remnant of Z3 symmetry in the 
pure gauge limit. In QCD with dynamical quarks, there appear 
three Z3 vacua, when T is larger than a critical temperature 
Te- The Z3 vacua are classified by the phase <p of <P, and 
each has anyone of 4>, <fi + 2ir/3 and <f> + 4tt/3. Roberge and 
Weiss [55] found that there is a first-order phase transition at 
9q = 7r/3 mod 27r/3 where the ground state is changed from a 
vacuum to its Z3 images; the RW phase transition is illustrated 
in Fig.|6](shown later). The transition is called the "Roberge- 
Weiss transition." In this transition, charge conjugation (C) 
symmetry is spontaneously broken and # q -odd quantities such 
as the phase ipof^ are order parameters of the transition i50Tl . 

In the ordinary PNJL model with the scalar-type four-quark 
interaction only, the chiral transition occurs at higher T than 
the deconfinement transition, unlike LQCD data at zero and 
imaginary /i q . In Ref. 114911 . we revealed that the PNJL model 
with the scalar-type eight-quark interaction 1 29, 30, 47], 



(16) 



G sS [(qqY + (qi-y 5 TqY 

and the vector-type four-quark interaction lT28i[30ll48ll . 

G v (q lv q)\ (17) 

can reproduce LQCD data at imaginary [i q . Since the cou- 
pling constants, G s % and G v , of the interactions are adjusted 
to the LQCD data, the correlation between a and <P is still 
weaker in this model than in LQCD. We then propose another 
possibility to explain the strong correlation shown in LQCD. 
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The origin of the four-quark vertex G s is the one-gluon ex- 
change diagram between two quarks and its higher-order dia- 
grams. If the gluon field A v has a vacuum expectation value 
(Aq) in its time component, A v is coupled to (Ao) which is 
related to <P through (O [57]; see Fig.Q]for the diagrammatic 
description. Hence, G s is changed into an effective vertex 
G S (<P) that can depend on <P J13l- The effective vertex G s (<£) 
is called the entanglement vertex, and all interactions includ- 
ing G s (<£) are referred to as the entanglement interactions. It 
is expected that the <P dependence of G s (^>) will be deter- 
mined in the future by an exact method such as ERG ll56l - 
IHt]. In this paper, however, we simply assume the follo wing 
G s that preserves chiral symmetry, C symmetry ll32l l50ll 
and extended Z3 symmetry 114711 : 

G S (<Z>) = G s [l - ai<£<r - a 2 (<£ 3 + <2>* 3 )]. (18) 

In the mean field approximation, the mesonic potential Um is 
modified as follows, 

U M (a,$) = G s [l -a^<P* -a 2 (<£ 3 +<P* 3 )]a 2 (19) 

and the constituent quark mass is changed into 

M = m -2G s (<2>)cr. (20) 

Thus, this model has entanglement interactions between a and 
<P in addition to the covariant derivative in the original PNJL 
model. The gap equation (0 can be evaluated by using the 
chain rules even in the presence of the entanglement interac- 
tions. 




Fig. 1: The diagrammatic description of the effective vertex G s ($). 

In this paper, the original PNJL model is simply called 
PNJL. The PNJL model with the entanglement vertex G s (<£) 
is referred to as entanglement PNJL (EPNJL), while the PNJL 
model with the scalar-type eight-quark and the vector-type 
four-quark interaction is referred to as PNJL-8V. 



such as the coincidence and the tuq dependence are preserved 
in the parameter region a\,ai ~ 0.20 ± 0.05. The validity 
of the parameter set in the EPNJL model is confirmed for real 
and imaginary fx\. 



A. Transitions at zero and finite quark-number chemical 
potentials 

First, we consider the case of /i q = 0. Figure [2] shows the 
T dependence of the order parameters a and <P, while Fig. [3] 
represents chiral and Polyakov-loop susceptibilities, Xa and 
X$, as a function of T, where the susceptibilities are normal- 
ized by T to become dimensionless 041 l53ll . As shown in 
Fig. |2j the chiral and deconfinement transitions are crossover 
in both the PNJL and EPNJL models. Figure^a) presents x<j 
and x<p m the PNJL model. The peak position of Xa, i-e., the 
critical temperature T a of the chiral transition, is much larger 
than the peak position of x<z> , that is, the critical temperature 
T$ of the deconfinement transition. Figure Ob) corresponds 
to Xa and x<s> m the EPNJL model. In this model, the two 
transitions coincide with each other within numerical errors. 
Thus, the entanglement vertex G s (<?) makes the correlation 
between the chiral restoration and the deconfinement transi- 
tion stronger, as expected. 
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Fig. 2: (color online). T dependence of the chiral condensate and the 
Polyakov loop at (9 q = 0. The curves that decrease (increase) as T 
increases represent the chiral condensate (Polyakov loop). The solid 
(dashed) curves are the results of the EPNJL (PNJL) model. Here, 
the chiral condensate is normalized by the value ctq at T = 0. 



III. NUMERICAL RESULTS 

In this section, we consider the case of Nf — 2, and take 
mo = 5.5 MeV unless otherwise mentioned. In the PNJL and 
PNJL-8V models, we take the same parameter set as in the 
previous analysis of Ref. [49]. In the EPNJL model, we take 
the same parameter set as the PNJL model, but To is taken 
to be 190 MeV so as to reproduce LQCD data at y« q = 0, 
when(ai,a2) = (0.2,0.2). This parameter set reproduces 
LQCD data at zero /j, q on the coincidence between T a and 
T$ IH |6jl |62|] and at imaginary on the m dependence 
of the order of the RW endpoint IU9I1 . Qualitative properties 



Next, we consider the case of 9 q — tt/3. Figure Hlpresents 
the T dependence of a and the absolute value of <P. In the 
PNJL model, the deconfinement transition at T = 189 MeV 
is first order, while the chiral transition is crossover; a has a 
small jump at T — 189 MeV, but it is just a discontinuity in- 
duced by the first-order deconfinement transition in \<P\. In the 
EPNJL model, the deconfinement transition at T = 185 MeV 
seems to be very weak first-order, since \<P\ has a small jump 
there within the present numerical accuracy, although it is not 
explicitly seen in Fig. |4] 

Figure |3a) shows that T a > T$ in the PNJL model, while 
Fig-Hb) shows that T a w T<g in the EPNJL model. Thus, the 
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Fig. 3: (color online). T dependence of the susceptibilities of the 
chiral condensate (dashed curve) and the Polyakov loop (solid curve) 
at 6 q = 0. Panels (a) and (b) correspond to the PNJL and EPNJL 
models, respectively. 



Fig. 5: (color online). T dependence of the susceptibilities of the 
chiral condensate (dashed curve) and the Polyakov loop (solid curve) 
at 9 q — 7r/3. Panels (a) and (b) correspond to the results of the PNJL 
and EPNJL models, respectively. 



entanglement vertex yields a stronger correlation between the 
chiral and deconfinement transitions also at 6 q = 7r/3. 
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Fig. 4: (color online). T dependence of the chiral condensate and the 
Polyakov loop at 8 q — n/3. The meaning of the curves is the same 
as in Fig. [2] 

Figure [6] shows the phase diagram in the 6 q -T plane. In 
the original PNJL model, T a is much higher than T$, while 
both are close to each other in the EPNJL model. The vertical 
dot-dashed lines at 9 q = ir/3 mod 2ir/3 are the RW transition 
line and the Z3 images. The endpoint of the RW transition 



line is located at T = T E ~ 189 MeV in the PNJL model and 
at 185 MeV in the EPNJL model. On the RW transition line 
at T > Te, C symmetry is spontaneously broken. As a con- 
sequence of this fact, 6> q -odd quantities such as the phase ip 
of the modified Polyakov loop 'P are discontinuous, while 8 q - 
even quantities have a cusp there B47I450I1 . Thus, the 6> q -odd 
quantities are order parameters of the RW phase transition. In 
the original PNJL model, the transition is second order ifsoh 
for the Polyakov potential proposed by Fukushima l33ll . but 
first order 15 ll 15211 for the Polyakov potential proposed by 
RoBner, Ratti and Weise [38]. The latter is more consistent 
with LQCD data at imaginary /i q than the former [51]. In the 
latter, the deconfinement phase transition is first order near the 
RW endpoint; the endpoint of the first-order deconfinement 
transition line is second order, and susceptibilities of several 
quantities diverge simultaneously there [51]. In the EPNJL 
model, such a first-order deconfinement transition line does 
not appear or very short if it does emerge, since the deconfine- 
ment transition at the RW endpoint seems to be a very weak 
first-order transition, as mentioned above. 

Figure [7] shows results of the EPNJL model for the RW 
phase transition. Panel (a) presents the T dependence of the 
phase il> of & at 8 q = ir/3 for the three cases tuq = 5, 150, 
and 400 MeV. The RW transition at the endpoint is first or- 
der for too = 5 and 400 MeV, but second order for too — 
150 MeV. In the limit of large toq, the transition is obviously 
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Fig. 6: (color online). Phase diagram in 8 q -T plane. Panel (a) is 
the result of the standard PNJL model with no entanglement vertex, 
while panel (b) is the result of the EPNJL model with («i, Q2) = 
(0.2,0.2). The solid (dashed) curves represent the deconfinement 
(chiral) transition. The vertical dot-dashed lines denote the RW tran- 
sition lines. Lattice data are taken from Ref. fl4ll . 
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Fig. 7: (color online). The RW phase transition in the EPNJL model. 
In panel (a), the phase of the modified Polyakov loop at t9 q = ir/3 is 
plotted as a function of T for three cases of light, intermediate, and 
heavy quark masses. Panel (b) shows the phase diagram of the RW 
phase transition in the mo-T plane. The solid (dotted) curve shows 
that the RW phase transition on the boundary is first (second) order. 



-0.0375 < nl < 0.08 [GeV 2 ] is expressed as 



first order, since the quark contribution to fi is suppressed 
and hence the deconfinement transition is controlled by the 
Polyakov potential U. Meanwhile, the RW endpoint is al- 
ways first order in the original PNJL model [51] and in the 
PNJL-8V model. Panel (b) shows the phase diagram of the 
RW phase transition in the mo-T plane; C symmetry is spon- 
taneously broken above the curve, while it is preserved below 
the curve. The solid (dashed) curve shows that the RW phase 
transition is first (second) order on the boundary. The critical 
mass mo(l — > 2) [mo (2 — > 1)] from the first-order (second- 
order) to the second-order (first-order) transition is rather sen- 
sitive to the numerical accuracy. In the present numerical ac- 
curacy, the critical masses are mrj(l — >• 2) = 50 ± 5 MeV and 
mo (2 — > 1) = 180 ±5 MeV. This mo dependence of the order 
of the RW endpoint is consistent with the recent result lfl9lr20ll 
ofLQCD. 

Figure [8|a) shows the phase diagram in the whole fi^-T 
plane obtained by the EPNJL model. The solid, dotted, dashed 
and dot-dashed curves represent the first-order chiral phase 
transition, the crossover chiral transition, the crossover de- 
confinement transition and the RW transition, respectively. 
The transition line (dashed and solid curves) in the region 



T = c + cx/j, 



(21) 



where c = 0.173 [GeV], c x = -0.377 [GeV" 1 ], and 
c 2 = -2.71 [GeV -3 ]. Point E is an endpoint of the RW tran- 
sition, while point C is a CEP of the first-order chiral phase 
transition. Point A is a meeting point between the RW tran- 
sition line and the crossover chiral transition line, while point 
B stands for the critical temperature of the chiral and decon- 
finement transitions at zero /u q . Locations of these points are 
tabulated in Table Ull Thus, there exists a CEP not only in the 
PNJL-8V model J49[] but also in the EPNJL model. 



A 13 C E 

(0-/3 x 193, 193) ( , 173) (160, 161) (in/3 x 185, 185) 

TABLE II: Locations (/i q , T) of points A, B, C and E. All locations 
are shown in MeV. 



Figure[8jb) shows the first-order chiral phase transition line 
and its CEP in the original PNJL, the PNJL-8V, and the EP- 
NJL models. The locations of the CEP in the three models 
are summarized in Table [Till The CEP is located at smaller 
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/i q and larger T in the EPNJL model compared with the other 
models. Thus, the entanglement vertex yields a drastic effect 
on the phase diagram at real fi q . 
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0.4 



Fig. 8: (color online), (a) Phase diagram in the fi q -T plane in the 
EPNJL model. The left (right) half-plane corresponds to imaginary 
(real) /i q . See the text for definitions of lines and points. Lattice data 
are taken from Ref. ffll , (b) The first-order chiral phase transition 
line and its CEP in the original PNJL, the PNJL-8V, and the EPNJL 
models. 



TABLE III: Summary of locations (fiq, T) of CEP in three models. 
All locations are shown in MeV. 



B. Transitions at finite isospin chemical potential 

The parameter set in the EPNJL model was determined 
in the previous subsection so as to reproduce LQCD data at 
zero and imaginary fi q . The validity of the parameter set is 
confirmed in this subsection for real and imaginary fi\ where 
LQCD data are available. 

The quark-number and isospin chemical potential, /i q and 
Mi, used in this paper are defined by 

flu + fJ-d flB flu — Md Miso 

Mq = —j— = T , Ml = — 5— = — (22) 



with the ?!-quark (d-quark) number chemical potential fi u 
(fid)- Here, fiQ an d Miso ®ce the baryon and original isospin 
chemical potentials coupled, respectively, to the baryon 
charge B and to the isospin charge I3. For comparison with 
LQCD, we use ^1 as the isospin chemical potential instead of 
the original definition Miso- 

The formalism of the PNJL model at finite mi is straight- 
forward from Sec. [TT] The only essential difference is that the 
pseudoscalar condensate ir = (qi^Tiq) is nonzero, in gen- 
eral, at finite Mi- Therefore, the E ± (p) in (|6]l are replaced 
by 



E±(p) = ^{E(p)+fi l f+N 2 ± 



for the u quark and 

E ± (p) = v/(£(p)-w) 2 +^ 2 ±Mq 



(23) 



(24) 



for the d quark, where N = —2G s ('P)ir. The mesonic poten- 
tial Um is a lso changed into 



U M = G a ($)(a 2 + ir 2 ). 



(25) 



PNJL 


PNJL-8V 


EPNJL 




(327, 124) 


(261, 124) 


(160, 161) 





See Refs. f5lll and ll53ll for the details of the formalism with 
finite mi! the only difference from the formalism is that G s is 
replaced by G s (<£) in the EPNJL model. 

First, we consider imaginary mi = idiT. In Fig. |9l we 
show T dependence of a and <P at 6\ = ir/2. In the standard 
PNJL model with no entanglement vertex, the critical temper- 
ature T a of the crossover chiral transition is about twice the 
critical temperature T,p of the first-order deconfinement tran- 
sition. This weak entanglement between the chiral restoration 
and the deconfinement transition still persists in the PNJL-8V 
model also [51]. 

The origin of the weak entanglement is the following. The 
u-quark loop contribution to the T-dependent part of Q is 
nearly canceled by the rf-quark loop contribution. Because 
of this cancellation, the thermal part i? th of Q is reduced at 
61 = 7t/2 and /i q = to 



-AT 



d :i 



p r 

(2tt) 3 I 



Tr n 



ln(l + ie-^ (p) -**) 
ln(l + } +U 



AT |i!P J {Tr c ln(l + L 2 e- 2 ^P))}+//. 



(26) 



because L — L' at // q = 0. In the last line of (|26l l, the ex- 
ponent in the first term is not f3E(p) but 2f3E(p), indicating 
that the temperature effect is reduced effectively by 1/2 in the 
first term. In the case of no entanglement vertex, the T depen- 
dence of a is controlled by the first term, while that of <P is 
controlled by U. Therefore, T„ w 2T$ in the original PNJL 
and the PNJL-8V model. 

In the EPNJL model, the entanglement vertex appears not 
only in the first term of the last line of (|26T l but also in the 
vacuum part of Q [the first and the Um terms in ©]. This 
induces a strong correlation between the chiral restoration and 



8 



PNJL 



PNJL-8V 



EPNJL 
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Fig. 9: (color online). T dependence of the chiral condensate and the 
Polyakov loop at 6\ — -k/2 and /i q = 0. See Fig.f2]for the meaning 
of lines. 



the deconfinement transition. Actually, as shown by the solid 
curves in Fig.|9j both the transitions are first order and T a = 
T<2>. LQCD data at 9\ — n/ 2 are not available in the two-flavor 
case but in the eight-flavor case [24]. The result of the EPNJL 
model is consistent with the LQCD result. 
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Fig. 10: (color online), (a) Phase diagram in the /j,i-T plane at /x q = 
in the EPNJL model. See the text for definitions of lines. LQCD 
data are taken from Ref. [23]. (b) The locations of the TCP at /i q = 
in the original PNJL, the PNJL-8V, and the EPNJL models. 

Next, we consider real fij. FigurefTUl'a) shows the phase di- 
agram in the fii-T plane at /i q = 0. The solid and dotted lines 



(401,171) (320,169) 



(95,170) 



TABLE IV: Summary of locations (jj,i,T) of TCP at fi q 
models. All locations are shown in MeV. 



in three 



stand for the first-order and second-order pion-superfluidity 
transitions, respectively. The meeting point between the solid 
and dotted lines is a TCP by definition. The crossover chiral 
and deconfinement transitions agree with each other, as shown 
by the dashed line. The EPNJL result reproduces LQCD re- 
sults on the chiral and deconfinement transitions and also on 
the pion-superfluidity transition. FigurefTOl'b) shows locations 
of the TCP at /z q = in the original PNJL, the PNJL-8V, and 
the EPNJL models. The entanglement vertex G S (<P) largely 
affects the location of the TCP. The locations of the TCP in 
the three models are summarized in Table II VI 



IV. SUMMARY 

In summary, we have extended the PNJL model by intro- 
ducing an entanglement vertex G s (<£) phenomenologically. 
The effective vertex generates entanglement interactions be- 
tween a and <P. The EPNJL model with G s (<£) can repro- 
duce two phenomena simultaneously; one is the strong cor- 
relation between the chiral restoration and the deconfinement 
transition that appears in LQCD at imaginary /i q and real and 
imaginary fij, and the other is the quark-mass dependence of 
the order of the RW endpoint predicted by LQCD very re- 
cently Ill9[|20ll , Thus, the EPNJL model is consistent with all 
LQCD data at imaginary /i q and real and imaginary /xt. 

The functional form of the entanglement vertex G s (<£) is 
determined by respecting extended Z3 symmetry, chiral sym- 
metry and charge conjugation symmetry. The strength of the 
entanglement vertex is determined by LQCD data at imagi- 
nary /i q , and the validity of this model building is confirmed 
by LQCD data at real and imaginary /zt. The entanglement 
vertex largely changes the location of the TCP in the /ii-T 
plane and the location of the CEP in the /ii-T plane. 

The present phenomenological approach seems to be com- 
plementary to the exact renormalization-group approach. It is 
highly expected that the functional form and the strength of 
the entanglement vertex will be determined in the future by 
the theoretical approach. 
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